arXiv: 1507.04547v2 [astro-ph.HE] 17 Jul 2015 


Mon. Not. R. Astron. Soc. 000, [THTT] (2015) Printed 20 July 2015 (MN style die v2.2) 


Attenuation of super-soft X-ray sources by circumstellar 
material 


M.T.B. Nielsen^*, M. Gilfanov^’^ 

^ Max-Planck Institut fiir Astrophysik, Karl-Schwarzschild-Str. 1, Postfach 1317, D-85741 Garching, Germany 
^ Space Research Institute of Russian Academy of Sciences, Profsoyuznaya 84/32,117997 Moscow, Russia 


Accepted -. Received 20 July 2015; in original form - 


ABSTRACT 

Recent studies have suggested the possibility of significantly obscuring super-soft X- 
ray sources in relatively modest amounts of local matter lost from the binaries them¬ 
selves. If correct, then this would have explained the paucity of observed super-soft 
X-ray sources and would have significance for the search for single-degenerate type 
la supernova progenitors. We point out that earlier studies of circumbinary obscura¬ 
tion ignored photo-ionisations of the gas by the emission from the super-soft X-ray 
source. We revisit the problem using a full, self-consistent calculation of the ionisa¬ 
tion state of the circumbinary material photo-ionised by the radiation of the central 
source. Our results show that the circumstellar mass-loss rates required for obcuration 
of super-soft X-ray sources is about an order of magnitude larger than those reported 
in earlier studies, for comparable model parameters. While this does not entrirely 
rule out the possibility of circumstellar material obscuring super-soft X-ray sources, 
it makes it unlikely that this effect alone can account for the majority of the missing 
super-soft X-ray sources. We discuss the observational appearance of hypothetical ob¬ 
scured nuclear burning white dwarfs and show that they have signatures making them 
distinct from photo-ionised nebulae around super-soft X-ray sources imbedded in the 
low density ISM. 

Key words: binaries: close - supernovae: general - white dwarfs - X-rays: binaries 
- radiative transfer - circumstellar matter - winds, outflows 


1 INTRODUCTION 

The topic of a massive white dwarf accreting hydrogen- 
rich material from a non-degenerate companion and pro¬ 
cessing that material into heavier elements on its surface 
is of relevance to the understanding of a number of as- 
troph ysical objects, including symbiotics (e.g. iMikolaiews 
I 2 OI 2 I) and single-degenera te type la supernova progeni¬ 
tors (e.g. iMaoz et al.ll2014l l. In the the case of the single¬ 
degenerate scenario, the accretion process is necessary for 
the white dwarf to grow in mass until it becomes mas¬ 
sive enough to undergo thermonuclear runaway in its de¬ 
generate core and explode. The steady burning of ac¬ 
creted hydrogen-rich material on the surface of a mas¬ 
sive white dwarf is believed to be poss ible only in a 
narrow interval of m ass-accretion rates llNomotd 1 19821 : 
IShen fc Bildsteiill2007l ). and the process is expected to emit 


raw 


of magnitude between the expected and observed number 
of systems, as well as the integrated luminosities of super- 
soft X-ray sources in old galactic populations, provided 
type la supernova progenitors are in fact si ngle-degenerates 
llGilfanov &: Bogdhnll2010l : lOi Stefandl201(]|) . 

The missing super-soft X-ray sources could perhaps be 
explained if the X-ray emitting systems were obscured by 
local material. There is growing evidence for the presence of 
circumstellar material around at least a subse t of type la su¬ 
pernova, 85 inferred from supernova spectr a (|Gerardy_eLalJ 
200d; I Patatjd^lj JgOOTl : ISinion et al.l 120091 : ISternberg et ahl 
2 OI 1 I: ICao et all I 2 OI 5 I) , as well as from historical super¬ 


nova remnants ([ Borkow ski et al.ll200d : IChiotellis et al.ll2012l : 
iBroersen et al ]l201^ 

Circumstellar material around type la supernova pro¬ 
genitors could be supplied in a variety of ways, e.g. by a 


However, several of the companion, incomplete accretion of wind-Roche lobe 

one to two orders overflow ijMohamed & Podsiadlowskil 2007 

) or tidal interac- 

tions between the components 1 Chen et al 
of circumstellar obscuration on the obser 

I 2 OIID. The effect 
vability of super- 
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soft X-ray emission from progenitor systems was investi- 
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gated by iNielsen et alj ll2013t) : IWheeler fc Poole'^ ll2013ll . 
Both studies suggested that relatively modest amounts 
of circumstellar material could obscure super-soft X-ray 
sources to the point where they would no longer be de¬ 
tectable with current X-ray telescopes. However, both stud¬ 
ies either ignored or used rather simple assumptions concern¬ 
ing the effect of photo-ionisation of the circumstellar mate¬ 
rial, and specifically did not take photo-ionisations of met¬ 
als into account. The present study seeks to remedy those 
shortcomings by conducting a full treatment of the photo¬ 
ionisation of the obscuring circumstellar material, to get a 
more accurate handle on how much circumstellar material 
is required for obscuration. 

The reader should note the terminology used in the fol¬ 
lowing: mass-transfer rate denotes the rate of mass trans¬ 
ferred from the donor to the accretor, regardless of whether 
the transferred material is retained on the accretor, whereas 
mass-accretion rate - or just accretion rate - is the rate 
of material transferred and retained on the accretor; lastly, 
mass-loss rate refers to the rate of mass lost from the bi¬ 
nary into the circumbinary region where it may obscure the 
super-soft X-ray source. We use the symbol for the lat¬ 
ter. 

In Section [2] we describe our model. Section [3] presents 
the results of our calculations, while Section 3] discusses the 
implications of the results. Section [5] concludes. 


2 MODEL 

The geometry of the model e xamined in th i s stud y is iden¬ 
tical to the one considered in INielsen et al.l (l2013l l. The pri¬ 
mary difference between that study and the present one is 
the incorporation of full photo-ionisation calculations. In 
particular, we consider a super-soft X-ray source inside a 
spherically symmetric shell of material, surrounded by vac¬ 
uum. Throughout the following, we assume that the source 
is a blackbody, with effective temperature fcTeff between 30 
eV and 150 eV, and luminosities Lboi between 10®^ erg/s and 
10 ^® erg/s, parameters charact eristic of ’canonical’, persis¬ 
tent super-soft X-rav sources |van den Heuvel et alJ 1 19921 : 
iKahabka fc van den Heuvellll997l. 

The gas shell surrounding the source is defined by an 
inner and outer boundary (rinn and rout, respectively), and 
the density is given by an r~^-law, characteristic of spheri¬ 
cally symmetric, constant-velocity wind mass loss: 

/ \ ^ ^ /"IN 

pv) = -.—5-> rinn < r < rout (1) 

where p(r) is the mass density as a function of radius, 
is the mass-loss rate into the circumstellar region and 
is the velocity of the wind, assumed to be constant in the 
region of interest. Inside the inner and outside the outer 
boundaries we assume vacuum. The inner radius can be en¬ 
visioned as equivalent to the orbital radius of the binary (see 
also Sections[3lfc l4.2ll . Realistic density profiles are likely to 
be more complex than this, and clearly depend on the mass- 
loss history of the binary, however, the simple r~^ profile is 
a useful first approximation and also facilitates comparison 
with previous results. 


2.1 Order-of-magnitude estimates 


The total mass of the gas shell is 

Mu, = ^ (runt-rinn) (2) 

Uuj 

= 9.5 X 10^® ( 1 - — ) M ™,-6 Mmo r-ont.Au g 

\ ^out J 

where is the mass loss rate in units of 10 ~® M 0 yr“^, 

uioo is the outflow velocity in units of 100 km/s and rout.Au 
is the outer radius of the shell in astronomical units. 

The column density of the gas shell is: 


Nh 


Mw 

47rMn,Mrn,p 


1 


rinn 


1 

^out 


( 3 ) 


1.4 X 10^^ 



'^100 


AU 


where /r « 1.4 for Solar abundance, and rrip is the mass of 
the proton. A column density of ~ 10^^ cm“^ of neutral ma¬ 
terial of Solar abundance can easily obscure emission from a 
su per-soft sourc e with t emperature of ~ 10 0 eV, a s discussed 
in [Nielsen et al.l (l2013l l: lwheeler fc PoolevI (l2013l l. However, 
photo-ionisation of the wind material by the emission from 
the super-soft source can significantly reduce its opacity. 

Similar to the classical Stromgren sphere problem, the 
effect of the photo-ionisation can be estimated, to the 
first approximation, by comparing the total rates of photo¬ 
ionisation and recombination in the gas shell: 


Ni = 


Nr = 


L„dv = f -— 

h,.o ^’^0 


n CKrec d r — 


Ml 


47™2,/i^mp 


1 

rinn 


1 

r out 


( 4 ) 

( 5 ) 


where hvo is the photo-ionisation threshold of the main ab¬ 
sorbing ion, OnEc its recombination coefficient, L is the lumi¬ 
nosity of the central source and the factor / depends on its 
radiation spectrum. In Eq. ®, it was assumed that the gas 
shell is optically thick to the radiation of the central source 
(in the absence of photo-ionisation). If the photo-ionisation 
rate Ni exceeds the recombination rate Nr, the shell material 
will be fully ionised and therefore transparent to the radia¬ 
tion of the central source. The critical value of the mass out¬ 
flow rate is determined by the condition Ni = Nr- At lower 
Mro the gas shell is fully ionised and transparent, whereas 
at higher outflow rates it is opaque to the radiation of the 
central source. For monoatomic hydrogen gas, the value of 
the critical outflow rate can be computed straighforwardly; 
with huo = 13.6 eV and Orec = 2.6 x 10“^® cm®sec“^ (case 
B) one obtains: 

Mu,,crit = 2.2 X 10-®(/L38)^/^uioo 

X fl-—Mo/yr (6) 

\ ^OUt J 


where Las is the luminosity in units of 10®® erg/s. For the 
central source temperature of fcTeff = 100 eV, for example, 
/ « 5 X 10“® and the leading factor in Eq. becomes 
4.9 X 10“® M 0 yr“^. Eq. ([ 6 ll is sufficiently accurate for hy¬ 
drogen gas shell, but gives only an order of magnitude es¬ 
timate for gas of solar abundance. The latter case requires 
a more careful account of the ionisation state of the main 
absorbing species. However, Eq. predicts a nearly correct 
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scaling of the critical mass loss rate (as confirmed by Cloudy 
calculations in the following sections): 

M™.crit oc r\lf (l - (7) 

\ out / 

As one can see from Eq. 0 , the critical mass loss rate rela¬ 
tively weakly depends on the luminosity of the central source 
and the inner radius of the shell and is nearly insensitive to 
the outer radius of the shell (for sufficiently large rout/rinn). 
The latter should be expected, given the steepness of the 
assumed density profile. 

The above considerations strictly apply to a stationary 
gas shell. In the case of the wind, (presumably) neutral ma¬ 
terial is continuously added at the inner edge of the shell. It 
is promptljQ photo-ionised by the radiation from the central 
source. However, if the new material is added at a sufficiently 
large rate - i.e. one comparable to the photon injection rate 
- it will screen the rest of the shell from the ionising radia¬ 
tion, thus reducing its degree of ionisation and increasing its 
opacity. The importance of advection is determined by the 
ratio of the rate at which new material is added to the inner 
edge of the shell to the injection rate of ionising photons: 

= 6 X if-iLss)-^ (8) 

^phot 

where /_i is the factor / in units of 10~^. The ratio given 
in Eq. m becomes of the order of unity at ~ 10 ^ 
M 0 yr“^, which is much larger than the mass outflow rates 
considered here. 

2.2 Cloudy calculations 

For detailed photo-ionisation calculations, we used the one¬ 
dimensional photo-ionisation code Cloudy, version 13.02 
(for a detailed descript ion of the latest release of Cloudy 
see iFerland et al.l l2013f ). Cloudy calculates the ionisation 
and thermal state of a cloud of gas that is irradiated by 
a central source, by simultaneously solving equations of 
photo-ionisation and thermal equilibrium. It is assumed that 
atomic processes in the gas have become time-steady. This 
assumption is reasonable, at least in the context of single¬ 
degenerate type la supernova progenitors, since all micro¬ 
physical time scales in the shell are much shorter than time 
scales relevant to super-soft sources. In our simulations, we 
fix the density distribution and compute the gas tempera¬ 
ture self-consistently, as maintained by the photo-ionisation 
of the gas by the central source. For the gas density distribu¬ 
tion, we use the built-in dlaw wind option in Cloudy which 
reproduces the density structure in Eq. ©, and the sphere 
keyword as appropriate for a closed geometry. The advec¬ 
tion effects are not included in our calculations, as justified 
above. Cloudy’s output files specify the emergent spectrum, 
temperature, pressure, and ionisation structure, in addition 
to a host of other data. 

We run Cloudy for three values of u™: 10, 50, 100 km/s. 
This spans the likely wind velocities expected from red giant 
and asymptotic giant star donors llOwoc 5 l2013li . For the 
purposes of this study we do not specify the nature of the 

^ the photo-ionisation time at the inner edge of the gas shell is 
of the order Tpf^ ~ 10“® — 10“^ s. 


wind. For each combination of Lboi, fcTefj, rinn, font and 
Uw we run Cloudy for a range of values of from 10“^° 
MQyr”*^ to 10 “® MQyr~'^. 

The chemical composition of t he ma terial in the shell is 
Solar, as given in ICrevesse et ahl (l2010l i. In S ubsection 13.31 
we co mpare this with the abundances used in iNielsen et al.l 
(|2013h . 

For hydrogen and helium, Cloudy can be expected to 
give reliable results up to n ~ 10^® cm“®. However, for num¬ 
ber densities 10 ^° cm“® the full treatment of a plasma 
containing heavier elements should not be considered reli¬ 
able, see the Cloudy manual ’Hazy0. For comparison, the 
largest value of for which one of our plots have data 
points (10“® MQ/yr at Vinn = 1 AU and Uw = 50 km/s on 
Figure |9]) corresponds to the number density of 1.9 x 10^° 
cm“®, i.e. within the nominal density range of Cloudy. 


3 RESULTS 

Cloudy’s output can be used to reconstruct the spectrum 
of the system consisting of the super-soft X-ray source and 
circumstellar cloud, along with the ionisation structure of 
the circumstellar material. In the following, we use the term 
output spectrum to mean the sum of the spectrum of trans¬ 
mitted radiation and the emission from the irradiated cloud. 
For the output data, we concentrate on the photon energy 
interval between 0.3 and 1.5 keV. Below 0.3 keV, even a mod¬ 
est ISM column effectively blocks any super-soft emission, 
and above 1.5 keV the emitted energy of a super-soft X-ray 
source is negligible. 

To quantify how the attenuation depends on the var¬ 
ious input parameters, we calculated the integrated emit¬ 
ted super-soft luminosity L of the output spectrum for each 
input configuration and compared this with the integrated 
emitted luminosity of the unabsorbed super-soft source spec¬ 
trum, Lo (using the photon energy interval between 0.3 and 
1.5 keV as integration limits). This yielded a transmission 
ratio L/Lo between 1 (no obscuration) and 0 (complete ob¬ 
scuration) for each input configuration. 

3.1 Attenuation of super-soft emission 

To demonstrate the general trends in our results, Figure [T] 
shows a comparison between the unabsorbed input spectrum 
of a fcTefi = 50 eV, Lboi = 10^® erg/s blackbody source, and 
the output spectra in a cloud with inner and outer radii of 
1 and 10 AU, respectively, and w™ = 50 km/s, for three 
values of from 1 x 10 “^ M 0 yr“^ to 1 x 10 “® M 0 yr“^. 
For the smallest mass-loss rate shown the output spectrum is 
only weakly attenuated. For larger values of M^, absorption 
edges of various neutral and ionised species appear in the 
spectrum along with their line and recombination emission. 
Among the absorption edges, the most significant is Ovil 
around -739 eV. At M = 1 x 10“® Mgyr"^ the 0.3-1.5 keV 
flux is attenuated by a factor of 0.30. 

The behaviour shown in Figure^is qualitatively similar 

^ Sect.3.8 on p.35 in ’Hazy’ Part 2 for Cloudy cl3.1. The two 
Hazy manuals are downloadable with the Cloudy installation files 
on http://nublado.org/wiki/DownloadLinks 
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Figure 1. Attenuated spectra of the super-soft X-ray source (blackbody, kT^g = 50 eV, -Lbol = 10®® erg/s) in the X-ray band for three 
values of the mass-loss rate. Parameters of the gas shell are rjun = 1 AU, rout = 10 AU, Um = 50 km/s. The dashed line shows the 
unabsorbed spectrum. Compare with the T.0-10 AU’ curve in the second panel from the top, middle plot in Figure|2] 


for other configurations of kT^g, rinn, Tout and u-w, although 
the critical mass-loss rate where obscuration becomes sig¬ 
nificant is of course different for different combinations of 
parameters. Results of our calculations are summarised in 
Figure [2] where we show the attenuation curves as a function 
of the mass loss rate for different values of the source tem¬ 
perature and gas shell parameters. Table [T] gives the M-w,ait 
values required to produce an attenuation of a factor of 10^. 
It also lists the corresponding column densities of the gas 
shell Nh- 

3.2 Parameter dependence 

Figure |3] shows comparisons between attenuation curves for 
different combinations of inner and outer radii, as a func¬ 
tion of the circumstellar mass loss rate. In agreement with 
the scaling suggested by Eq. ©, the attenuation curve is 
moderately dependent on the value of the inner radius and 
nearly insensitive to the position of the outer boundary of 
the shell, as long as rout/rinn is sufficiently large. 

Dependence of the attenuation curve on the luminosity 
of the central source is illustrated in Figure [4] As one can 
see from the plot, it is somewhat weaker than predicted by 
Eq. 0. This is because the derivation of Eq. © did not 
take collisional ionisation into account. 

As the gas density in the shell is determined by the 
ratio M^juw, the attenuation curves shift along the 
axis linearly with the outflow velocity. This is obvious from 
Figure [2] and Table [T](cf. Eq. 0). 



Figure 3. Attenuation curves showing the effect of different val¬ 
ues of rout (top panel) and rinn (bottom panel). In both panels, 
fcTefj = 50 eV, Lbol = 10®® erg/s, = 50 km/s. 


3.3 Chemical abundances 

The chemi cal abundances used in Nielsen et al. I iI2oT 3) were 
taken from I Anders fc E biharal ll 19821) ( for all other elem ents 
than helium) and Peimbert fc Torre s-P eimbert (Il977) (fo r 
helium), as tabulated in iMorrisoii fc McCammon 1 19831) . 
The present study used the more update d chemical abun¬ 
dances table from iGrevesse et al. I ll2010l) . Figure [5] shows 
a comparison between the attenuation curves for the two 


abundance tables, when calculated with Cloudy using other¬ 
wise identical input parameters. For reference, the ’GASSIO’ 
curve in Figure [5] is the same as the ’1.0-10 AU’ curve in 
the second panel from the top, middle column on Figure [2] 
As Figure shows, the difference between the transmission 
curves computed for the two abundance tables is insignif¬ 
icant. For comparison, a transmission ratio of 10“'* corre- 
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Figure 2. Attenuation of a super-soft source as a function of the mass outflow rate Mw’ The super-soft radiation is assumed to have 
a blackbody spectrum with luminosity of Lboi = 10^^ erg/s, and the luminosity emerging from the system is integrated from 300 eV 
to 1500 eV. The three columns of panels show results for different outflow velocities: Uw = 10 km/s (left) Uw = 50 km/s (middle) 
and Uw = 100 km/s (right). Different panels in each column correspond to different values of the source temperature (marked in the 
upper-right corners of each panel) and different curves in each panel correspond to different values of the inner radius of the gas shell, 
as marked in the plots. The outer boundary of the shell is rout = 10 AU in all cases. 




Figure 4. Attenuation curves for different values of the lumi¬ 
nosity of the central source. Other parameters are fixed at the 
following values:/cTeff = 50 eV, Uw = 50 km/s. 

sponds to a mass-loss rate of 2.04 x 10“^ M 0 yr“^ for the 
MM83 abundances and ~ 2.18 x 10“^ M 0 yr“^ for GASSIO. 


Figure 5. Comparison b etween attenuation cu rves computed for 
the abundance table of l Grevess^^^a]jj2^^ (marked GASSIO) 
and the older one used in NielseiTet al.l ll2013h (marked MM83), for 
/cTeff = 50 eV, Lboi = 10^® erg/s, rinn = 1 AU, rout = 10 AU 
and Uw = 50 km/s. 


4 DISCUSSION 

The results presented above demonstrate the importance 
of the proper account for photo-ionisation when consider¬ 


ing attenuation of super-soft X-ray sources by circumbinary 
material. In particular, the values of Mw required for signif¬ 
icant obscuration of a canonical super-soft X-ray source are 
about an ord e r of m agnitude larger than those reported in 
iNielsen et al.l (|2013l) . Figure [ 6 ] shows a direct comparison be- 
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Table 1. Values of the mass-loss rate and corresponding Nu, 
needed for an attenuation of -Lobs/^O = 0.01. 


Uyj 

Teff 

^’inri) ^’out 

^'W,CY\t 

Nh 

[km/s] 

[eV| 

[AU| 

[M0yr-l| 

[cm“2] 

10.0 

30.0 

0.1, 10 

3.6x10“® 

5.2x10^2 

II 


0.5, 10 

1.3x10“’^ 

3.2x1022 

II 

” 

1.0, 10 

2.0x10“’^ 

2.6x1022 

II 

50.0 

0.1, 10 

8.2x10“® 

1.2x1023 

II 


0.5, 10 

2.1x10“’^ 

5.8x1022 

II 


1.0, 10 

3.5x10“'^ 

4.5x1022 

II 

100.0 

0.1, 10 

1.6x10“’^ 

2.3x102® 

II 

” 

0.5, 10 

3.5x10“'^ 

9.5x1022 

II 


1.0, 10 

5.1x10“’^ 

6.6x1022 

II 

150.0 

0.1, 10 

2.0x10“’^ 

2.8x102® 

II 

” 

0.5, 10 

4.1x10“’^ 

1.2x1023 

II 

” 

1.0, 10 

6.1x10“’^ 

7.9x1022 

50.0 

30.0 

0.1, 10 

1.8x10“’^ 

5.2x1022 

II 

” 

0.5, 10 

5.9x10“’^ 

3.2x1022 

II 

” 

1.0, 10 

1.0x10“® 

2.6x1022 

II 

50.0 

0.1, 10 

4.1x10“’^ 

1.2x1023 

II 


0.5, 10 

1.1x10“® 

5.8x1022 

II 

” 

1.0, 10 

1.7x10“® 

4.5x1022 

II 

100.0 

0.1, 10 

8.0x10“’^ 

2.3x1023 

II 


0.5, 10 

1.8x10“® 

9.5x1022 

II 

” 

1.0, 10 

2.6x10“® 

6.6x1022 

II 

150.0 

0.1, 10 

9.9x10“’^ 

2.8x1023 

II 

” 

0.5, 10 

2.2x10“® 

1.2x1023 

II 


1.0, 10 

3.1x10“® 

7.9x1022 

100.0 

30.0 

0.1, 10 

3.7x10“'^ 

5.2x1022 

II 


0.5, 10 

1.2x10“® 

3.2x1022 

II 

” 

1.0, 10 

2.0x10“® 

2.6x1022 

II 

50.0 

0.1, 10 

8.2x10“’^ 

1.2x1023 

II 


0.5, 10 

2.2x10“® 

5.8x1022 

II 


1.0, 10 

3.5x10“® 

4.5x1022 

II 

100.0 

0.1, 10 

1.6x10“® 

2.3x1023 

II 

” 

0.5, 10 

3.5x10“® 

9.5x1022 

II 

” 

1.0, 10 

5.1x10“® 

6.6x1022 

II 

150.0 

0.1, 10 

2.0x10“® 

2.8x1023 

II 


0.5, 10 

4.3x10“® 

1.2x1023 

II 

” 

1.0, 10 

6.1x10“® 

7.9x1022 


tween the two works. In the figure, we plot the attenuation 
of the integrated number of photons emitted by a n obscured 
source as given in Table 1 of iNielsen et al.l ll2013ll and those 
found in the current study for the same input parameters. As 
already explained above, the difference in these two curves 
stems from the detailed treatment of photo-ionisation of the 
circumbinary material by the radiation from a super-soft 
source. 

The significantly larger mass-loss rates required for ob¬ 
scuration of the super-soft source provide a stronger con¬ 
straint on the possibility of nuclear burning massive white 
dwarfs being hidden by absorption by l ocal circumstellar 
material. For example, in iNielsen et al.l ll2013ll it was sug¬ 
gested that for a 1 and 10 AU inner and outer radius, and a 
wind velocity of 10 km/s it would be possible for a massive, 
steady-burning white dwarf accreting at a few times 10 ~^ 
M 0 yr“^ to be more or less fully obscured if < 10 % of the 
material lost from the donor star did not become accreted 
onto the white dwarf but instead were available for obscura¬ 
tion in the circumstellar region. In the current study, we find 



Figure 6. Compar ison between attenu ation of number of ionising 
photons (N/Nq) in iNielsen et al.l ll2013ll and the current study, for 
fcTeff = 50 eV, Lbol = 10^* erg/s, rinn = 1 AU, rout = 10 AU, 
and Unj = 10 km/s. 

for the same input values of parameters, that this amount of 
obscuration would require a mass loss into the circumstellar 
region of > 50% of the mass-transfer rate, i.e. the amount of 
material lost into the circumstellar region would have to be 
of the order as the material being accreted onto the white 
dwarf. For wind velocities of ~ 50 — 100 km/s the circumstel¬ 
lar mass loss rate should exceed > 90% of the mass-transfer 
rate, i.e. ~ 10 times more material should be lost into the 
circumstellar space than be accreted on to the white dwarf. 
Recalling that the steady-burning rate for a massive white 
dwarf is a few times 10 ~^ M 0 yr“^, this means that to sus¬ 
tain both the steady burning super-soft X-ray source and 
the obscuring envelope the donor would have to lose > 10 ~® 
M 0 yr“^ for a > 10® years. The high mass-loss rate and the 
total mass lost (around a Solar mass or more) constrains the 
population of possible donors, and hence binary systems ca¬ 
pable of producing significantly obscured, super-soft X-ray 
sources that eventually reach the Chandrasekhar mass and 
becomes type la supernovae. 

Although we considered here a specific case of the cir¬ 
cumbinary wind, interpreted broadly, our results show that 
a canonical super-soft source is capable of ionising about 
Nh ~ 10 ^^ — 10 ^® cm~^ of material located at the distance 
of a ~ few astronomical units from the source (Table [T]). 

Looking from a different perspective, a narrow but 
dense shell of gas with the hydrogen column density in ex¬ 
cess of ~ 10 ^® cm~^ around a single-degenerate type la su¬ 
pernova progenitor could conceivably obscure the super-soft 
X-ray emissions from such a system. Because of the large ve¬ 
locities of supernova ejecta (upwards of ~ 10,000 km/s - or 
5.8 AU/day - for type la supernovae) such a shell would be 
overtaken by the ejecta within a matter of days. Hence, up¬ 
per limit measurements from e.g. radio observations would 
need to be made almost immediately after the supernova 
explosion to provide useful constraints on the amount of 
circumbinary material present at the ~AU distances from 
the supernova. The example of nearby SN2011fe is proba¬ 
bly to be considered a best case in this respect, as it was 
discovered and subjected to follow-up observatio ns almost 
immediately after exploding llNueent et al.ll201lli . Even so, 
the radio limits placed on SN2011fe were based on observa- 
tions conducted a litt le over two days after the supernova 
llChomiuk et al.l I2OI2I ’) , which means that the blast wave 
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Figure 7. Broad band (infrared to X-ray) spectra of an attenuated super-soft X-ray source for three different attenuation levels; an 
unobscured source (M = 10~® MQ/yr), moderately obscured source (M = 5 X 10“^ M 0 /yr) and heavily obscured source (M = 3 X 10~® 
Mo/yr). The model parameters are at their default values: fcTgff = 50 eV, Lbol = 10^® erg/s, rinn,rout =1 AU,10 AU. The dashed black 
line is the incident spectrum. 


would have moved to a distance of ~ 12 AU (assuming the 
above-mentioned ejecta velocity) at the time observations 
began. Hence, a very dense shell of material of, say, between 
1 and 10 AU from the exploding object may conceivably 
have eluded detection, as the supernova ejecta would have 
passed through the obscuring medium before observations 
began. More distant type la supernovae are typically dis¬ 
covered several days to weeks after the explosion, and in 
such cases radio observations should not be expected to yield 
upper limits capable of ruling out narrow configurations of 
circumstellar material. However, since single-degenerate pro¬ 
genitors of type la supernovae need to be super-soft X-ray 
sources for extended periods of time and require continuous 
mass-loss rates > 10“®Myr~^ to remain hidden, the extent 
and mass of the circumstellar gas surrounding such sources 
could be much larger than the narrow configurations men¬ 
tioned above and should therefore, theoretically at least, be 
detectable via post-supernova radio observations. Also, such 
extended systems should be detectable in pre-supernova ob¬ 
servations in the recombination li nes of hydrogen and helium 
and forbidden lines of metals (cf. IWoods fc Gilfanov|[2013l ). 


4.1 Observational appearance of obscured SSS 

Although circumbinary obscuration is unlikely to explain the 
paucity of super-soft sources in the context of progenitors 
of type la supernovae, it is plausible that stably nuclear- 
burning white dwarfs experience various degree of circumbi¬ 
nary obscuration. In particular, this may be one of the fac¬ 
tors (along with the reasons intrinsic to the population syn¬ 
thesis calculations) explaining the discrepancy between the 
number of super-soft sources predicted in pop ulation syn¬ 
thesi s stud ies and their observed numbers (e.g. IChen et al.l 
|2014 I 2 OI 5 I ). While quenching the bulk of the X-ray emis¬ 
sion from the stably nuclear-burning white dwarf, circumbi¬ 
nary obscuration leads to the appearance of various emission 
signatures at all wavelengths (including X-ray band itself). 


which can be exploited to search and identify such obscured 
super-soft X-ray sources. 

Figure [7| shows the output spectrum of an attenuated 
super-soft X-ray source in the broad wavelength range from 
infrared to X-ray wavelengths. Spectra for three values of 
the mass-loss rate are shown, corresponding to differ¬ 
ent obscuration levels in the X-ray band - from completely 
unobscured to completely obscured source (cf. Fig[2]). In the 
completely unobscured case, the gas in the wind is nearly 
fully ionised and optically thin, and the attenuation at these 
mass-loss rates is negligible (save for a short period when 
the source initially ’turns on’ and ionises the cloud). Such 
sources would be similar to naked super-soft sources when 
observed in X-rays, as discussed previously. In the obscured 
case, the soft X-ray radiation is reprocessed to the optical 
and UV bands, giving rise to copious emission lines. In the 
intermediate case of moderate obscuration a number of emis¬ 
sion features can be observable in the X-ray band as well. 


4 . 1.1 X-ray lines 

X-ray lines, although not the brightest among the emis¬ 
sion features in our calculations, can carry away up to 
— 10“® of the luminosity of the central source. The 
brightest among them are resonant lines of H-like ions and 
He-like triplets of the low-Z metals. In Figure 0 we show lu¬ 
minosities of several of the brightest X-ray lines as a function 
of the mass-loss rate. As the amount of the circumbinary ma¬ 
terial increases with the mass-loss rate, luminosities of these 
lines initially grow with Mw, until obscuration in the X-ray 
region sets in at ~ 10“® M©/yr (see Figure[2|. At this point, 
the X-ray line luminosities plummet due to absorption by 
the neutral outer layers of the wind. The behaviour of par¬ 
ticular lines is mainly determined by the changes of relative 
abundances of respective ions, with the intrinsic line ratios 
being characteristic for photo-ionised plasma with ionisation 
parameter of ^ « 37.4 Las MZ^ uioo- However, as the mass- 
loss rate approaches the critical value, the escaping line lu- 
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Figure 8. Luminosity of the most prominent X-ray emission lines vs the mass-loss rate, for our default model setup {kT^g = 50 eV, 
^bol = 10®® erg/s, Uw = 50 km/s, rinn,rout =1 AU,10 AU). 


minosities and their ratios are significantly modified by the 
absorption by the neutral outer part of the wind, giving rise 
to the patterns shown in FigO 

iNess et ^ l|2ni,'j l studied high-resolution spectra for a 
number of super-soft X-ray sources in local galaxies. Based 
on these observations, they suggested a subdivision of super- 
soft X-ray sources into two classes: SSa (’absorbed’) and 
SSe (’emitting’) - the former having generally blackbody¬ 
like spectra with some prominent absorption features, and 
the latter showing strong emission features (comparable in 
strength to the continuum) superimposed on weak black- 
body spectra. When correlated with the inclinations of the 
systems of that study (and with the caveat that incli¬ 
nation determination is quite problematic in most of the 
cases) members of the SSa class were predominantly high- 
inclination systems, while SSe systems were predominantly 
viewed more edge-on. This led Ness and collaborators to 
suggest that the spectral differences of the two classes might 
arise as a result of viewing-angle, rather than from intrin¬ 
sic differences between the systems themselves, conceptually 
similar to the viewing-angle effect in the unified model for 
AGNs. Although the direct comparison with observed spec¬ 
tra is outside the scope of the present work, our study does 
recreate several o f the features of the SSe class reported in 
iNess et"^ I ll2013ll in the interval where the central sources 
transition from unobscured to fully obscured regime. The 
two strongest emission features in the sp ectra of U Sco and 
CAL 87 on Figure 3 in iNess et al.l (l2013l l correspond to the 
O Vlll-doublet at 18.97A, 18.98A, and the line blend con¬ 
sisting of the N Vll-doublet at 24.78A, 24.79A, and the N 
VI at 24.89A. We hnd the same features are prominent in 
our spectra, as seen on our Figure [T] Several other, less pro¬ 
nounced lines in our spectra can be also found in the ob¬ 
served spectra of U Sco and CAL 87, e.g. Ne IX at 13.45A, 
O VIII at 21.6A, and possibly three N VI lines around 29A 
(28.78A, 29. lA, 29. SA). At lower X-ray energies th e simi¬ 
larities between our spectra and those of iNess et al.l ll2013l l 
become less obvious, presumably due to interstellar absorp¬ 
tion in the observed spectra. The fact that our relatively 


simple model setup can rep roduce som e of th e features of 
the SSe spectra disucssed in iNess et al.l (l2013l) is not in it¬ 
self surprising, since any model with significant incident lu¬ 
minosity in X-rays which interacts with a Solar metallicity 
gas is likely to produce a somewhat similar set of emission 
features. However the luminosities of these lines, their ratios 
and the presence and strength of absorption lines and edges 
are determined by the details of the model such as geom¬ 
etry of the gas and emission pattern, ionisation parameter 
and chemical abundances. Detailed analysis of these effects 
is beyond the scope of the present paper. 


4-1-2 Optical lines 

As the mass-loss rate grows towards full obscuration, 
the X-ray emission is gradually quenched, and the absorbed 
energy increasingly re-emitted at UV and optical energies. 
In this regime, copious emission lines in the UV, optical 
and infrared bands are produced in the photo-ionised wind 
material. The majority of these lines (except for extreme UV 
lines), unlike X-ray lines, escape the wind material without 
significant attenuation and therefore may serve as a valuable 
diagnostic tool in the search for and identification absorbed 
super-soft sources. 

Figure shows the luminosity of a number of obser- 
vationally relevant lines of hydrogen and helium as well 
as forbidden lines of metals as a function of the mass-loss 
rate. The luminosity of hydrogen and helium recombination 
lines grows monotonically with mass-loss rate and then sat¬ 
urate as the central source becomes completely obscured (at 
roughly ~ few x 10 “® M 0 /yr in our default model, 
see Figure [21). At this point, the ionised medium becomes 
confined within the neutral outer layers of the wind, and 
forbidden lines begin to appear as well, produced at the in¬ 
terface of the ionised and neutral gas. At Mu, > 3 — 5 x 10“® 
M©/yr, the Lja luminosity of the gas cloud reaches ~ 10% 
of t he total luminosit y of the central source. 

iRappaport et all lll994h modelled the effect of per¬ 
sistent, unobscured, super-soft X-ray sources on the sur- 
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Figure 9. Luminosity of several most prominent optical and UV 
lines vs the mass-loss rate for our default model (fcTeff = 50 eV, 
Lboi = 10^^ erg/s, Uw = 50 km/s, rinn^'^out =1 AU,10 AU). The 
’wavy’ behaviour of the Ha and curves at Mw ^ 10 ~® M 0 /yr 
is a numerical effect, arising from optical depth convergence prob¬ 
lems in Cloudy at large optical depths. This is a known problem 
that will presumably be addressed in an upcoming update of the 
code (Cloudy Collaboration, private communication). 


Figure 10 . Comparison between the [OIII] 5007A and Hell 
4686A to H^ line ratios expected from the ionised nebulae in the 
low density ISM (Rappaport et al. (1994), red triangles) and the 
corresponding ratios expected for an obscured super-soft source 
imbedded in the dense wind considered in this work (blue circles). 
The wind case was computed for our default model (fcTeff = 50 
eV, L = 10^^ erg/s, rinn??"out = 1-0,10.0 AU); the points shown 
are for Myj values for which the [OIII] luminosity is non-negligible, 
i.e. Mw ^ 3.5 X 10“® M 0 /yr. 


rounding (low density) interstellar medium. Their results 
showed that such sources should develop parsec-sized ioni¬ 
sation nebulae which should exhibit strong emission in for¬ 
bidden optical lines, most importantly [O III] 5008A, as 
well as in recombination lines of hydrogen and helium, in¬ 
cluding He II 4686A. They suggested using forbidden oxy¬ 
gen lines ([01] and [OIII]) as a diagnostic tool to find 
ionised ISM nebulae around on unobscured super-soft X - 
ray sources dRappaport et a0ll994l : iDi Stefano et al]|l995ll . 
Much higher gas densities in the wind 10® cni~® in the 
wind vs. ~ 1 cm~® in the ISM) make forbidden lines much 
weaker in our case. Figure [10] shows a comparison of the 
[OIII] 5 OO 7 A and Hell 4686A to line ratios in the case 
of the obscuration by the wind wi th those expected from 
the ionisation nebulae computed bv iRappaport et al.l (Il994l l 
(their Figure 5). As one can see from the plot, a super-soft 
source obscured by the dense wind considered in this work 
will be easily distinguishable from the case of the ionised 
nebula around a super-soft source imbedded in the low den¬ 
sity ISM0 

To conclude, the observational signatures of a fully ob¬ 
scured super-soft X-ray sources should be significant (up to 
10 ®^ erg/s for a 10^® erg/s source) Ly^ emission, cou¬ 
pled with the low [OIII]-to-H^ ratio, of ~ 10“®, and rela¬ 
tively high Hell 4686A-to-H;3 ratio, ~ 1. These character¬ 
istics make obscured super-soft sources an observationally 
distinct class of compact line emission sources. 

4.2 The inner wind radius 

The critical mass-loss rate at which obscuration of the cen¬ 
tral source becomes significant shows moderate dependence 

® A caveat here is that such a nebula has been so far de¬ 
tected only around one super-soft source llRemillard et aLlll995t 
IWoods &: Gilfanovll201^ . 


on the assumed inner radius of the wind. Eq. [6] predicts 
that Mw.crit oc whereas numerical calculations with 

Cloudy show marginally stronger scaling. The moderate 
fimi-dependence is a consequence of the two assumptions 
made in our calculations: (i) that the pgaa oc density pro¬ 
file extends to infinitely small radii and (ii) that the source of 
radiation is located in the locus of this density distribution, 
i.e. at the point r = 0. These two assumptions combined 
together lead to the strong dependence of the total recom¬ 
bination rate in the gas shell on its inner radius, eq. ®. 
To eliminate these assumptions, one would need to specify 
the driving mechanisms behind the mass loss and accurately 
consider the physical conditions at the inner wind boundary, 
which is far beyond the scope of this work. Below we con¬ 
sider the example of symbiotic systems to attribute some 
physical meaning to the inner radius in our si mulations. 

Recent modelling of the red giant winds lISuzu ^ l2007f) 
demonstrated that the pgas oc density profile may in¬ 
deed be a valid representation of the wind structure down 
to rather small radii, as low as ~0.1 AU. The source of ionis¬ 
ing radiation (i.e. the white dwarf), however, is offset by the 
distance a (binary separation) from the singular point of the 
gas distribution. In the absence of the singularity, the depen¬ 
dence of the total recombination rate in the gas shell on its 
inner radius is insignificant (as long as rinn rout) and the 
only role of the rinn in this case is to provide the normal¬ 
isation of the density distribution, via the mass continuity 
equation Mu, = A-k p{r) u-w ■ Therefore it is natural to set 
rinn = a in this case, i.e. the inner radius is of the order of 
the orbital separation in the binary system. Note also, that 
the 3-dimensional nature of the problem becomes essential 
at r ~ a, as the red giant will be casting a shadow onto the 
wind material located behind it. Obviously, the simplified 
1-dimensional consideration in this paper is applicable only 
to the lines of sight which are not passing close to the donor 
star. 
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To conclude, while the picture of a spherical shell con¬ 
fined between the inner and outer boundaries is a simplifi¬ 
cation of the real geometry of a circumbinary gas cloud, our 
model is a useful first approxim ation. This is the g eometry 
considered in the previous work (iNielsen et al.ll2C)13 1 and its 
use here facilitates comparison with those results and helps 
to isolate the effect of photo-ionisation. Future work (cur¬ 
rently in progress) will aim at dealing with more realistic 
density structures and source spectra than those discussed 
here. 


5 CONCLUSIONS 

We have presented full photo-ionisation simulations of the 
interaction between the radiation from a super-soft X-ray 
source and a circumstellar cloud of Solar metallicity gas. 
This seeks to model a single-degenerate system that has 
lost material from either the donor or the accretor into the 
circumstellar region. The aim has been to determine the 
amount of material (parametrised by a circumstellar mass- 
loss rate) required to significantly obscure the X-ray emis¬ 
sion from the accretor. 

To simulate the effect of the X-ray source on the gas 
cloud w e used the one-dim ensional photo-ionisation code 
Cloudy (iFerland et al.|[2013h . Our analysis explored inter¬ 
vals of source temperature, luminosity, inner and outer radii 
of the gas cloud, and wind velocity, to cover the range of 
values for realistic systems. 

We found that the amount of circumstellar matter re¬ 
quired to hide a super-soft X-ray source from our current 
observational capabilities is about an order o f magnitude 
larger than that found in INielsen et all (l2013ll . In particu¬ 
lar, the required mass-loss rates from the binary typically 
exceeds 10”^ to a few times 10~® Meyr”^ for the parame¬ 
ters considered in the ’default’ case in this work. Since the 
donor star has to provide both the material being accreted 
onto the white dwarf (and hence powering its X-ray emis¬ 
sion) and the circumstellar material, this places significant 
constraints on the mass of the donor star. 

We discussed the observational appearance of hypothet¬ 
ical nuclear burning white dwarfs obscured by a wind. Such 
systems would appear as luminous (~ 10®^ erg/s) compact 
sources of Lya emission with strong Hell recombination lines 
(Hell 4686A to H^ ratio 1 and suppressed forbidden 
lines of metals, e.g. [OHI]-to-H ,3 ratio, of 10“^), mak¬ 
ing them distinct from photo-ionised nebulae predicted to 
exist around super-soft X-ray sources imbedded in the low 
density ISM. 
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